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ABSTRACT 

The amplitude of weak lensing should increase with source distance, rising steeply 
behind a lens and saturating at high redshift, providing a model-independent means 
of measuring cosmic geometry. We measure the amplitude of weak lensing with redshift 
for three massive clusters, A370 {z = 0.375), ZwC10024+17 (z = 0.395) and RXJ1347- 
11 {z = 0.451), using deep, three-colour Subaru imaging. We define the depth of 
lensed populations with reference to the COSMOS and GOODS fields, providing a 
consistency check of photo-z estimates over a wide range of redshift and magnitude. 
The predicted distance-redshift relation is followed well for the deepest dataset, A370, 
for a wide range of cosmologies, and is consistent with less accurate data for the 
other two clusters. Scaling this result to a new survey of 25 massive clusters should 
provide a useful cosmological constraint on w, complementing existing techniques, 
with distance measurements covering the untested redshift range, 1 < 2 < 5. 

Key words: cosmology: observations - gravitational lensing - galaxies: clusters: 
individual(Abcll 370) - galaxies: clusters: individual(ZwCl 0024.0-1-1652) - galaxies: 
clusters: individual(RX J1347.5-1145) 



1 INTRODUCTION 

Constraining cosmological parameters has been the focus 
of major surveys in the last decade, via precision cos- 
mic microwave background (CMB) temperature correlations 
(Spergel et al. 2007; Brown et al. 2009) and SN-Ia light 
curves (Riess et al. 1998; Perlmutter et al. 1999). A stan- 
dard, ACDM, cosmological model has been defined by this 
work, albeit at the price of accepting an accelerating expan- 
sion driven by a cosmological constant, and non-baryonic 
dark matter (DM) of an unknown nature dominating the 
mass density of the Universe. Measurements of the angular 
diameter distance of the CMB refer to 2: ~ 1100, and lumi- 
nosity distances are derived from SN-Ia in the range z < 1. 



* Based on data collected at Subaru Telescope and obtained from 
the SMOKA, which is operated by the Astronomy Data Center, 
National Astronomical Observatory of Japan. 



In principle, lensing can provide a complementary distance 
measurements in the range, z > 1, from the purely geomet- 
ric defiection of light, which increases with source distance 
behind a lens. 

For lensing clusters, the bend-angle of light scales lin- 
early with angular diameter distance ratio, dis/ds, the sepa- 
ration between the lens and the source, divided by distance 
to the source. This distance ratio has a characteristic geo- 
metric dependence on redshift, rising steeply behind the lens 
and then saturating at large source redshift (e.g.. Fig. 1 of 
Broadhurst, Taylor & Peacock 1995). This effect has been 
detected behind massive lensing clusters, where the sepa- 
ration in angle between multiple images of higher redshift 
sources is noticeably larger than for lower redshift sources. 
For example, for the well studied cluster SDSSJ1004+4112 
(z = 0.68), five images of a QSO a.t z — 1.734 are within 
an Einstein radius of about 6e ~ 7" (Inada et al. 2003; 
Oguri et al. 2004), whereas a more distant multiply lensed 
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galaxy behind this cluster at z = 3.332 is at a much larger 
Einstein radius of 6e ~ 16" (Sharon et al. 2005). Other 
sets of multiple images show this increasing angular scal- 
ing with source redshift, following well the expected general 
form of the redshift-distance relation in careful strong lens- 
ing analyses of deep Hubble data (Soucail, Kneib & Golse 
2004; Broadhurst et al. 2005; Zitrin et al. 2009b; Zitrin & 
Broadhurst 2009; Zitrin et al. 2009a). 

However, these studies are not able to distinguish be- 
tween the relatively subtle changes between cosmologies in 
the range of interest, due to the inherent insensitivity of the 
distance ratio di^/ds to the cosmological parameters. More- 
over, the bend-angle is particularly sensitive to the gradient 
of the mass profile, requiring many sets of multiple images 
in the strong regime to simultaneously solve for both the 
cosmological model and the mass distribution. Instead in 
practice, strong lens modeling usually adopts the standard 
cosmological relation in order to better derive the mass dis- 
tribution, with multiply lensed sources forced to lie on the 
lensing distance-redshift relation. This helps eliminate the 
otherwise considerable degeneracy in constraining the slope 
of the lensing mass profile of a galaxy cluster (Broadhurst 
et al. 2005; Zitrin et al. 2009b; Zitrin & Broadhurst 2009; 
Zitrin et al. 2009a). The hope of constraining the cosmog- 
raphy from strong lensing data is remote with current tools 
(Gilmore & Natarajan 2009). 

Weak lensing (WL), by contrast, offers a model inde- 
pendent way of constraining the cosmological parameters 
via the distance-redshift relation (Taylor et al. 2004). Image 
distortions and magnification depend on gradients of the de- 
flection field and in the WL limit, these are just proportional 
to dis/ds- The mass profile enters only in the stronger regime 
as a second order correction (Medezinski et al. 2007). Here 
we are concerned with the amplitude of cluster WL depen- 
dence on source redshift, for which no mass reconstruction 
is required when evaluating the distance-redshift relation. 
However, although this observed effect is independent of the 
mass distribution in the weak limit, the sensitivity to cos- 
mological parameters is still inherently very small. 

The analogous effect in the field has been explored in 
terms of the cosmic shear, to measure the general mass dis- 
tribution. Optimal formalisms have been developed which 
cross-correlate the foreground distribution of galaxies along 
a given line of sight with the distribution of background 
images (Wittman et al. 2001; Jain & Taylor 2003; Bacon 
et al. 2003; Taylor et al. 2004, 2007) with clear detections 
of large scale structure, including the COMBO-17 fields 
(Brown et al. 2003; Kitching et al. 2007) and the COSMOS 
field (Massey et al. 2007; Schrabback et al. 2010). To usefully 
derive cosmological parameters from general cosmic shear 
work, deep all-sky surveys have been proposed (e.g., LSST^, 
DES^ JDEM^ EUCLID''). 

Here we make use of detailed colour-colour (CC) in- 
formation for three intermediate redshift lensing clusters, 
A370 {z = 0.375), ZwC10024-Hl7 {z = 0.395), RXJ1347-11 
{z = 0.451) and define several samples of galaxies of differing 

^ http://www.lsst.org/lsst 

^ https : //www . darkenergy survey . org/ 

^ http://jdem.gsfc.nasa.gov 

* http : //sci . esa. int/science-e/www/object/index . cf m? 
f objectid=42266 



Table 1. The Cluster Sample: Redshift and Filter Information 



Cluster 


z 


Filters used^ 


Seeing 






(exp' time in sec) 


(arcsec) 


A370 


0.375 


Bj (7200), Rc (8340), 2'(14221) 


0.6 


ZwCI0024-|-17 


0.395 


Bj(3600), Rc(5280), 2' (1680) 


0.8 


RXJ1347-11 


0.451 


yi(1800), Rc(2880), 2'(4860) 


0.76 



^Detection band marked in bold. 

background depths with which to explore the dependence of 
WL distortion on source redshift. With only three bands 
we cannot reliably define photometric redshifts for a sizable 
proportion of objects, but instead, by reference to the red- 
shift disributions of the well studied deep field surveys, we 
may reliably define several galaxy populations of differing 
mean depths in the CC-space covered by the filters used for 
each cluster. 

Three colour selection has also been applied in a simi- 
lar context to simulations aimed at forecasting the capabil- 
ities of WL tomograghy (Jain, Connolly & Takada 2007; 
Medezinski et al. 2010). These simulations convincingly 
demonstrate that greater efficiency is likely by using lim- 
ited 3-band imaging for WL tomography, rather than by in- 
vesting greater imaging time in additional bands to improve 
photometric redshift precision. For our purposes too we show 
here that a judicious choice of non-orthogonal boundaries in 
CC-space allows the definition of several distinct redshift 
samples of differing mean depth, with relatively little over- 
lap in redshift. 

Here we rely on well studied deep field surveys to esti- 
mate the redshift distribution of these different background 
populations, using the wide-field COSMOS 30-band photo- 
metric redshift survey (Ilbert et al. 2009) and also the deeper 
GOODS-MUSIC survey which has wide multi-wavelength 
coverage in 14 bands (from the U band to the Spitzer 8 /im 
band) (Grazian et al. 2006). 

In § 2 we present the cluster observations and data re- 
duction and in § 3 we explain the selection of background 
galaxy samples and in § 4 we describe the WL analysis and 
outline the formalism. In § 5 we derive the WL amplitude 
and mean redshift information of the background samples, 
presenting our results regarding the lensing distance-redshift 
relation. We discuss the requirements for constraining the 
cosmological model with this method in 6. We summarize 
and conclude in § 7. 



2 SUBARU DATA REDUCTION 

We analyze deep images of three intermediate-redshift clus- 
ters, A370 ZwC10024-M7 and RXJ1347-11, observed with 
the wide-field camera Suprime-Cam (Miyazaki et al. 2002) 
in several optical bands, at the prime focus of the 8.3m 
Subaru telescope. The clusters are publicly available from 
the Subaru archive, SMOKA^. Subaru reduction software 
(SDFRED) developed by Yagi et al. (2002) is used for flat- 
fielding, instrumental distortion correction, differential re- 

^ http://smoka.nao.ac.jp 
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fraction, sky subtraction and stacking. Photometric catalogs 
are created using SExtractor (Bertin & Arnouts 1996). 
Since our work relies much on the colours of galaxies, we pre- 
fer using isophotal magnitudes. We use the Colorpro (Coe 
et al. 2006) program to detect in the _Rc-band and measure 
colours through matched isophotes in the other two bands. 
Astrometric correction is done with Scamp (Bertin 2006) us- 
ing reference objects in the NOMAD catalogue (Zacharias 
et al. 2004) and the SDSS-DR6 (Adelman-McCarthy et al. 
2008) where available. The observational details are listed 
in Table 1. 



3 SAMPLE SELECTION FROM THE 
COLOUR-COLOUR DIAGRAM 

For each cluster we use Subaru observations in three broad 
optical passbands and all observations are of good seeing, 
representing some of the highest quality imaging by Subaru 
in terms of depth, resolution, and colour coverage. We first 
describe how we separate the background galaxies from fore- 
ground and cluster galaxies, combining WL measurements 
and the distribution of objects in the CC plane and their 
clustering relative to the center of the cluster. We then ex- 
amine the redshift distribution of objects selected to lie in 
the background using the CC plane with reference to the 
COSMOS field where deep photometric redshifts are estab- 
lished to faint limits using 30 independent passbands cover- 
ing a very wide range of wavelength. 

Note that with only three bands, only a small fraction 
of the objects have well defined photometric redshifts, in 
the sense of having a single, narrow peak in their probability 
p{z), and therefore their use to separate among different red- 
shift populations is very limited. Using the BRz' CC space 
in this way with reference to the now well established red- 
shift surveys is arguably more reliable in separating galaxy 
populations of differing depths, in agreement with the sim- 
ulations of Jain et al. (2007). 

In our previous analysis of these data, Medezinski et al. 
(2010, hereafter MIO), we demonstrated how the cluster 
and foreground galaxies can be reliably identified and sep- 
arated from background galaxies in the CC diagram, using 
Bj,Rc, z bands (Fig. 1). We found in the field of A370 the 
prominent overdensity of galaxies centered on _Bj — Rc ~ 2 
and _Rc — z' ~ 0.8 denotes the red-sequence of cluster galax- 
ies (see Fig. 1, left-hand panel, where this overdensity is 
enclosed by dashed white line). This was also shown by the 
relatively small mean distance from cluster center of galax- 
ies in that region in CC space (Fig. 1 in MIO). The WL 
measurements for this population of objects is very close to 
zero, with a measured tangential distortion profile, grir), 
consistent with zero all the way out to the virial radius, as 
expected for cluster galaxies which are unlensed. The main 
central overdensity around Bj — Rc ~ 1 and _Rc — z' ~ 0.3 
consists of many foreground galaxies (see Fig. 1, left-hand 
panel, with an overdensity enclosed by solid white line). 
These galaxies show a very low level gr compared to the ref- 
erence background, marked in gray on Fig. 1 (right panel), 
and their surface density profile shows only modest central 
clustering (see MIO) indicating that most of these objects 
lie in the foreground of the cluster. When selecting back- 
ground galaxies we stay well away from these regions of CC 



space which are dominated by cluster and foreground galax- 
ies to minimize contamination by these unlensed galaxies, 
as described below. 

To identify background populations in MIO, we took 
into account both the WL signal and the density distribu- 
tion of galaxies in the CC plane. In CC space a relatively red 
population can be rather well defined, dominated by an obvi- 
ous overdensity (around Bj — Rc ~ 0.5 and -Rc — z' ~ 0.8), 
and the bluest population is confined to a separate cloud 
(around Bj — Rc ~ 0.3 and _Rc — z' ~ 0.2) of faint galax- 
ies with a clear lensing signal. Looking at the WL profiles 
of the red and blue samples, we see very similar behaviour, 
with a continuous rising signal toward the cluster center. 
Both cases show good agreement with each other (MIO). 
Combined together, they form our reference background 
sample, to which all the other samples we derive below will 
be normalized. 

The validity of our selection was also demonstrated in 
MIO by comparison with the spectral evolution of galax- 
ies calculated with the stellar synthesis code Galev^ (Ko- 
tuUa et al. 2009). Here we overlay the CC diagram with 
colour-tracks of galaxy models: E-type (exponentially de- 
clining SFR with Zq), so (gas-related SFR with Zq), Sa 
(gas-related SFR with 2.5^q), and Sd (constant SFR with 
0.2Zq). These evolutionary tracks originate in the low- 
redshift cloud we established as foreground in CC space and 
evolve to higher redshift passing through the cluster redshift 
and then to bluer colours as shown in Fig. 5 (right), and fi- 
nally end at the top-left corner of our CC space, dropping 
out of the i3,i-band at a redshift of z ~ 3.5 

In this paper we add additional samples of background 
galaxies. Two samples will consist of galaxies at redshifts not 
far beyond that of the cluster, which we term "orange" and 
"green". These are selected to lie on the upper-right of the 
CC space, corresponding reasonably to colour-tracks of E/SO 
galaxies which are predicted to show a bend in the CC plane, 
becoming bluer {Bj — Rc ~ 2 — 2.5 and Rc — z' ~ 1 — 1.5) 
toward higher redshift, matching well our observed distribu- 
tion. A third group comprises the red-cloud described above, 
centered on Bj — Rc ~ 0.5 and Rc — z' ~ 0.8, we call 
the "red" sample. This sample consists of an overdensity of 
background galaxies from the known "red" branch of the 
bimodality of field galaxies colours (Capak et al. 2007). A 
Forth sample consists of the prominent blue peak identified 
as background galaxies, lying at redshifts beyond the red 
galaxies, corresponding to the "blue" branch of the colour 
bimodality, where the colour-tracks turn redder at Bj — Rc 
as the UV-flux starts to drop out of the blue Bj band with 
increasing redshift, above a redshift of approximately, z ~ 2. 
Finally, we select galaxies from the top-left corner of the CC 
plane, corresponding to the predicted location of _Bj-dropout 
galaxies at an average redshift of z ~ 3.5. 

For each sample selected we examine the colour bound- 
aries of each sample as a function of the WL signal, to ensure 
we keep well clear of foreground or cluster members which 
otherwise dilute the WL signal of background galaxies. This 
approach has been investigated in our earlier work, where 
this problem was first identified as a major problem for WL 
work and rectified using appropriate CC-selection (Medezin- 



http : //www. galev . org/ 
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ski et al. 2007; Umetsu & Broadhurst 2008; Medezinski et al. 
2010; Umetsu et al. 2010). 

We apply this selection scheme for two more clusters, 
ZwC10024+17 and RXJ1347-11, though here the data is not 
as deep. Furthermore, the colour coverage is less wide in the 
case of RXJ1347-11 (only VRz') allowing the selection of the 
corresponding orange, green, red and blue background sam- 
ples, as above, but not sufficient for separation of dropout 
galaxies for a significant WL detection. 

The greater depth of the A370 imaging allows us to fur- 
ther divide the green, red and blue background samples into 
independent bright and faint subsamples for our WL mea- 
surements. For each of these samples, we now estimate their 
average WL signal, by comparing the tangential distortion 
measurements to that of the reference background sample, 
and derive what we call mean "gx amplitude" (see § 5.1). 
Subsequently, we also estimate the median redshift from the 
COSMOS photo-2 catalogue and its equivalent dis/ds (see 
§ 5.2). We summarize the properties of the selected samples 
in Table 2. 



4 WEAK LENSING MEASUREMENTS 

To make the WL catalogs, we use the IMCAT package devel- 
oped by N. Kaiser^ to perform object detection and shape 
measurements, following the formalism outlined in Kaiser, 
Squires & Broadhurst (1995, hereafter KSB). Our analysis 
pipeline is described in Umetsu et al. (2010). We have tested 
our shape measurement and object selection pipeline using 
STEP (Heymans et al. 2006) data of mock ground-based 
observations (see Umetsu et al. 2010, § 3.2). Full details of 
the methods are presented in Umetsu & Broadhurst (2008) , 
Umetsu et al. (2009) and Umetsu et al. (2010). 

The shape distortion of an object is described by the 
complex reduced-shear, g = gi + ig2, where the reduced- 
shear is defined as: 



ga = 7^/(1 - k). 



(1) 



The tangential component gr is used to obtain the az- 
imuthally averaged distortion due to lensing, and computed 
from the distortion coefficients gi,g2- 



gr = — (<?i cos 29 + g2 sin 26) , 



(2) 



where 6 is the position angle of an object with respect to 
the cluster centre, and the uncertainty in the gr measure- 
ment is (jt = ag/y/2 = o" in terms of the RMS error ag for 
the complex shear measurement. To improve the statistical 
significance of the distortion measurement, we calculate the 
weighted average of gr and its weighted error, as 



E» ^g,i 9T,i 



(3) 



(4) 



where the index i runs over all of the objects located within 
the n-th annulus with a median radius of 0„, and Ug^i is the 



inverse variance weight for i-th object. 



^ http : //www. if a.hawaii . edu/~kaiser/imcat 



where is the softening constant variance. We choose 
Q = 0.4, which is a typical value of the mean RMS ag over 
the background sample. We accurately combine the photom- 
etry with weak-lensing measurements of as many galaxies as 
possible, discarding objects below the seeing limit (given in 
table 1) plus two standard deviation of that value in the 
detection band, to remove stars and avoid unreliable shape 
measurements. 



4.1 Formalism: Relative Distortion Strength 

For a given source redshift Zs and a fixed lens redshift zi, 
the observable (complex) reduced gravitational shear g{zs) 
in the subcritical regime is expressed in terms of the gravi- 
tational shear 7 and the lens convergence k as (e.g., Seitz & 
Schneider 1997; Medezinski et al. 2007) 



g(zs) = ^{Zs){l - k[zs]) ^ = 7oo ^ l3''^''{Zs 



(5) 



where Koo and 700 are the lensing convergence and the grav- 
itational shear, respectively, calculated for a hypothetical 
source at Zs — >■ 00, and P{zs) is the lensing strength of 
a source at Zs relative to a source at Zs — ^ 00, I3{zs) = 
D{zs)/D{zs — >■ 00); D{zs) = dis/ds. Hence, the reduced 
shear averaged over the source redshift distribution is ex- 
pressed as 



(g) = 7oo 

where {/?''') is defined such that 



_ Jdz,N{z,)l3'={zs) 



JdzN{z 



(6) 



(7) 



with the redshift distribution N{zs). In the WL limit where 
l/too|, I7I00 <C 1, then 



(g) ^ (/3)7oc = {7). 



(8) 



Thus, the mean reduced shear is simply proportional to the 
mean lensing strength, {/?) oc (D). The next order approxi- 
mation is 



(7) 



(9) 



where fp = {/?^) /{/?)^ is a redshift-moment ratio of the order 
of unity (Seitz & Schneider 1997). 

Since the tangential distortion signal, ^t, is a function of 
cluster radius, we first decompose the tangential distortion 
profile of our background sample (B), defined above as our 
reference (see § 3), into the following form: 



9t,b{0) = asO 



(10) 



where the 
power-law 
distortion 
tice, we fit 
regime 
taneously 
&B. Next, 
power law 
ai to vary 



radial shape of grid) is assumed to be a single 
with a power index bs, and as represents the 
amplitude of our reference background. In prac- 
the outer profile of grid), excluding the nonlinear 
<1), to the power-law model, constraining simul- 
the distortion amplitude ob and the outer slope, 
for each of our other defined samples, we fit a 
with the same slope be, but allow the amplitude 
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Table 2. CC-sclcctcd Sample Properties 



1.1 n citfiT' 


TYlTll 


TTl n CTTl it"! iri O llTni'hc; 
lllOjtilll L Uv_i.l_^ 111111 LO 


N 


n 


p 














Oil t^lliili 


(.9 'J' -amplitude ratio) 


(PL) 




A370 


foreground 


18 < z' < 22 


1235 


1.3 


-0.04 


13/6 


0.33 




orange 


19 < z' < 23 


450 


0.5 


0.41 


1/6 


0.73 




green-bright 


20 < z < 22.5 


839 


0.9 


0.85 


11/6 


0.9 




green-faint 


22.5 < z' < 25 


1240 


1.3 


0.88 


6/7 


0.99 




red-bright 


22 < z' < 24 


5454 


5.6 


0.87 


17/8 


1.11 




red- faint 


24 < z' < 26 


6986 


7.1 


1.1 


27/8 


1.14 




blue-bright 


23 < z' < 24.6 


1679 


1.7 


0.95 


6/8 


1.79 




blue-faint 


24.6 < z' < 25.5 


2857 


2.9 


1.15 


23/8 


1.78 




drops 


24 < z' < 26.5 


1529 


1.6 


1.3 


13/8 


3.86 




background 


22 < z' < 26 


19362 


19.8 


1 


19/8 


1.31 


ZwC10024+17 


foreground 


18 < z' < 22 


734 


0.9 


0.08 


5/9 


0.29 




orange 


18 < z' < 23 


687 


0.9 


0.55 


11/9 


0.73 




green 


21 < z' < 25 


1582 


2 


0.83 


23/9 


0.98 




red 


21 < z' < 25.5 


7716 


9.7 


0.97 


8/9 


1.11 




blue 


23 < z' < 25 


2420 


3 


0.99 


13/9 


1.62 




drops 


25 < z' < 27 


625 


0.6 


0.93 


23/9 


3.71 




background 


21 < z' < 25.5 


10488 


13.1 


1 


7/8 


1.24 


RXJ1347-11 


foreground 


19 < z' < 24 


2452 


3.4 


-0.04 


14/8 


0.43 




orange 


19 < z' < 26 


473 


0.5 


0.76 


5/8 


0.76 




green 


21 < z' < 25 


1296 


1.1 


1.1 


5/8 


0.94 




red 


21 < z' < 26 


4880 


3.9 


0.92 


9/8 


1.11 




blue 


23.5 < z' < 26 


1509 


2.1 


1.18 


35/8 


1.95 




background 


21 < z' < 26 


7139 


6.3 


1 


4/6 


1.24 





0.5 1 1.5 -0.5 0.5 1 1 5 2 

R-Z R-Z 

Figure 1. Le/i.- Number density in Bj — Rq vs. Kq — z' CC space for A370. The four distinct density peaks are shown to be different 
galaxy populations - the reddest peak in the upper right corner of the plots (dashed white line) depicts the overdensity of cluster galaxies, 
whose colours lay on the red sequence; the middle peak lying blueward of the cluster comprises mainly foreground galaxies (solid white 
line); the two peaks in the bottom part (bluest in Bj — Rq) can be demonstrated to comprise of blue and red background galaxies. 
Together, the red+blue galaxies will serve as our "reference" background sample for WL purposes. Right: Bj — Rq vs. Rq — z' CC 
diagram, showing the distribution of galaxies in A370. Marked are the selected foreground sample (gray) and background samples: 
orange (orange), green (green), red (red), blue (blue) and dropout (magenta) background galaxies, selected to include galaxies lying away 
from the cluster and foreground regions. Overlaid are synthetic colour tracks including evolution, calculated with the Galev code for an 
elliptical, SO, Sa and Sd type models. 



QtAS) = are-"'' . (11) 

Therefore, if we calculate the lensing signal of i-th sample 
relative to the reference background (B), 

r, = gTA9)/9T,B{S) = <^i/aB. (12) 

© 0000 RAS, MNRAS 000, 000-000 



From equation (8), we obtain the following expression in the 
WL approximation (Kk)|, |(7)| ^ 1): 

r, = a,/aB ^ (/3)./(/3)b = {D)^/{D)b, (13) 
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where {)i {i = 1, 2, B) represents averaging over the red- 
shift distribution Ni{zs) of i-th galaxy sample. The relative 
distortion strength Ti can be regarded as a function of the 
discrete background sample i with the redshift distribution 
Ni{zs), which is observationally available and calibrated by 
deep, multi-band blank surveys such as the COSMOS sur- 
vey. For a given cosmological model, one can readily con- 
struct its theoretical prediction Ti (i = 1,2, ...,B) using a 
set of redshift distribution functions Ni{zs). The function 
Vi can be formally labeled by its mean redshift 

{zs)i = j dzsNi{zs)zs/ J dzsNi{zs) (14) 

which is independent of the cosmological model. 

To the next order of approximation, the distortion am- 
plitude ratio is written as (Appendix B of Medezinski et al. 
2007) 



A370 



Ti = 



{D)b 



{1 + ih.m^ - //3.b(/3)b) K^{e) + o(a^L)} (15) 



with /^., = (/3')i/(/3)? and /;3.b = (/3')b/{/3)|. The 
next order correction term is proportional to (f^.i{l3)i — 
/;3,b{/3)b)koo (^), which is much smaller than unity for the 
galaxy samples of our interest in the mildly nonlinear regime 
{6 > 1). We thus simply adopt equation (13) obtained in 
the WL approximation. This can be further justified by the 
fact that the slope parameter 6b is constrained by a least 
fit to the outer distortion profile. The mean weighted 
cluster radius (9) = '^iUg,i9i/ ^^Ug,i used for a fit is 
(fi) ~ 10 — 11 arcmin for our clusters, where the weak lensing 
approximation is valid. 



8 [arcmin] 

ZwC10024 



6 [arcmin] 

RXJ1347 



5 RESULTS 

5.1 Weak lensing profiles 

For each cluster, a reference background sample has been 
defined (explained above in § 3). The WL tangential distor- 
tion, gr, vs. distance from cluster center, 6, of each reference 
background is plotted in Fig. 2 (black crosses). Each refer- 
ence background sample is fitted by a power-law according 
to Eq. 10, but only in the WL regime, i.e., outside 9 > 1'. 
The fit is estimated using two ways - we fit the entire sam- 
ple dataset, weighting each galaxy gr.i by Ug^i (see § 4), 
and we also fit the binned gr profile, {gxifin)), weighting by 
the bin error, 1/ ar"^ {9n)- We find good consistency between 
the two fitting schemes. The goodness-of-fit values of the 
binned fit are displayed next to each power-law fit, and also 
in Table 2. As can be seen, a simple power-law serves as a 
reasonable fit in all cases. 

For each of our defined samples (foreground, orange, 
green, red, blue and dropouts) we again plot gx vs. radius 
in Figure 4 (gray, orange, green, red, blue and dropouts, 
respectively, top to bottom panels of each cluster). We fit 
each sample with the same power-law index, feB, given by 
its relevant reference background (the sample fit is shown 
as dashed black line with 1 — a confidence bounds, and the 
reference background fit is also shown as a dotted line) . The 
WL amplitude of each sample relative to the reference back- 
ground is given in each panel as F; = Oi/as (as defined by 
Eq. 12) and detailed in Table 2. We see that indeed, for 



-4 



Figure 2. Tangential distortion qt vs. distance from cluster cen- 
ter for A370 (top), ZwC10024+17 (middle) and RXJ1347-11 (bot- 
tom) reference background samples. Overlaid is the power-law fit 
(dashed black line) with 1 — cr confidence levels (shaded region) 
and the PL-fit is indicated. 



the foreground samples, the gT{9) profile agrees with zero 
throughout, and gives a relative WL amplitude of zero. 

As another consistency check, we plot the galaxy sur- 
face number density vs. radius for A370 samples in Figure 3. 
As can be seen, no clustering is observed toward the center 
for any of the samples, which demonstrate that there is no 
contamination by cluster members in the samples compris- 
ing only background members. The foreground sample (gray 
triangles) shows a modest increase in number density (factor 
of 2 increase from 6^ = 20 to 6^ = 2) compared to the cluster, 
despite the exclusion of the cluster early-type galaxies by 
colour. Bluer later-type cluster members are to be expected 
here given the redshift window sampled by reference to Fig- 
ure 6, which shows that the tail of the distribution reaches 
just beyond the redshift of A370 § 5.2). 

The data for A370 represents the best available data. 
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Figure 3. Galaxy surface number density vs. radius for A370 
foreground sample (gray triangles), orange (orange diamonds), 
green (green pentagrams), red (red squares), blue (blue circles) 
and dropout (magenta hexagrams) background samples. 

both in terms of the filter coverage - B, R, z' , and in terms 
of total exposure times in each band. For RXJ1347-11 we 
only have V, 7?, z' coverage, making it somewhat harder to 
set apart the different populations, and also dropping out 
between V and R is not as clear since the two filters are less 
well separated in wavelength, preventing us from selecting a 
dropout sample. For ZwC10024+17 we do have B, R, z' , but 
the z' band is much shallower. 

5.2 COSMOS photometric redshifts 

To estimate the respective depths of the different samples 
defined above from our Subaru photometry, we make use of 
the accurate photometric redshifts derived for the well stud- 
ied multi-band field survey, COSMOS (Capak et al. 2007). 
For COSMOS, photometric redshifts have been derived by 
Ilbert et al. (2009) using 30 bands in the UV to mid-IR. 
Since the COSMOS photometry does not cover the Subaru 
Rc band, we estimate iZc-band magnitudes for it. For this 
we use the HyperZ (Bolzonella, Miralles & Pello 2000) tem- 
plate fitting code to obtain the best-fitting spectral template 
for each galaxy, from which the Rc magnitude is derived 
with the transmission curve of the Subaru 7?c-band filter 
(see Umetsu et al. 2010). 

We then select samples by applying the same 
CC/magnitude limits as we did above for each of our clus- 
ters - A370, ZwC10024-|-17 and RXJ1347-11. This is shown 
in Fig. 5 (left) for the COSMOS catalogue and plotted in 
terms of the same CC plane as A370. The colour distribu- 
tion of COSMOS field galaxies seen in this B,R, z' CC-space 
is very similar to that of A370, displaying the same mor- 
phology, including red, blue and dropout populations, but 
without the density peak associated with the massive clus- 
ter A370. We also show how redshift varies in this CC-plane 
by calculating the mean photo-z redshift from COSMOS in 
fine bins over the CC plane (Fig. 5, right), with the sam- 
ples boundaries displayed as well. This demonstrates that 
the main overdensity in the CC plane near Bj — Rc ~ 1 and 
Rc — z ^ 0.3 has a mean redshift of around z < 0.5, which 
agrees with our estimation for A370 where we found very lit- 



Figure 6. Redshift distribution of all A370 samples: foreground 
(dashed gray), orange (dotted-dashed orange), green (green) 
bright (solid) and faint (dashed), red (red) bright (solid) and 
faint (dashed), blue (blue) bright (solid) and faint (dashed), and 
dropout (dotted-dashed magenta) using COSMOS photo-z's. 




redshift 



Figure 7. Redshift distribution of the dropout-selected sample 
using COSMOS photo-z's (light red) and using GOODS-MUSIC 
photo-z's (dark blue). The low-z peak seen more notably from 
the COSMOS photo-z's is most likely due to misclassified galaxy 
redshifts, supported by the smaller numbers found when using 
the GOODS-MUSIC photo-z catalogue. 

tie WL signal implying these object lie predominantly in the 
foreground of the cluster. We also see from this figure that 
the region where we picked "red" galaxies corresponds to 
z ~ 1 — 1.5, and the "blue" galaxies occupy a region of mean 
redshift around z ~ 2. Most notably, the top left corner of 
"dropout" galaxies corresponds to high-z with z > 3.5. We 
further plot the redshift distribution of all the samples in 
Fig. 6. We calculate the median redshift of each sample and 
summarized in Table 2. 

However, if we look at the distribution of COSMOS 
redshifts of the dropout sample (Fig. 7, red), we find it is 
somewhat double-peaked, with most galaxies lying around 
z ~ 3.5, but a significant fraction identified as having 
z ~ 0.4. Since we are certain most of the galaxies in this 
region are in fact high redshift dropout galaxies, justified 
by the apparent high WL distortions measured for these 
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Figure 4. vs. cluster radius for A370 bright samples (top left: foreground, orange, green, red & blue) and faint (top right: green, 
red, blue & dropouts), ZwC10024+17 (bottom left: foreground, orange, green, red, blue & dropouts) and RXJ1347-11 (bottom right: 
foreground, orange, green, red & blue), where the fixed power- law fit is overlaid (dashed blaek line) with 1 — a confidenee levels (shaded 
region). Also plotted is the power-law fit of the equivalent "reference" background sample (dotted curve). In each case the resulting 
normalized qt amplitude ratio, F, is denoted next to the profile. 



galaxies (see Fig. 4) and by reference to deep spectroscopic 
work in this rising plume of "dropout" galaxies (Steidel et al. 
1999), we may conclude the low redshifts assigned to some 
of these galaxies may be misclassified as low redshift early 
type/dusty galaxies in the photo- 2: catalogue of the COS- 
MOS field. This is not surprising, since at faint magnitudes 
the COSMOS photo-z have a relative high catastrophic fail- 
ure rate (Ilbert et al. 2009). 

We further examine this issue using the somewhat 
deeper GOODS-MUSIC catalogue (Grazian et al. 2006; San- 
tini et al. 2009), which has 15 bands, including high qual- 
ity ACS photometry (GOODS-S) and deep IRAC imaging 
which is very helpful in reducing the outlier rate a high-z. 
Using the z'-band selected GOODS-MUSIC photo-z, we de- 
rive a spectral classification for each galaxy using a library 
of ~ 200 PEGASE (Fioc & Rocca-Volmerange 1997) tem- 
plates very similar to that described in Grazian et al. (2006) , 
as included in the EAZY software (Brammer et al. 2009), 
and then calculate Subaru magnitudes for all the galaxies 
using the BPZ (Bem'tez 2000) code. 

By making the same CC selection, we plot the redshift 



distribution of the same dropout sample (Fig. 7, blue), but 
here practically no low-z peak is observed, and all galaxies 
in this region are estimated to have high redshifts, z > 2.5. 
This comparison between COSMOS and GOODS-MUSIC 
redshifts allows us to securely set a conservative lower B — R 
limit to avoid inclusion of real low-z objects. We can thus 
safely assume all objects identified as low- 2: in the COSMOS 
catalogue are largely mistakenly classified, a point also made 
in relation to this by Schrabback et al. (2010). Estimating 
the median redshift of the dropout sample gives z ~ 3.8, 
a value very close to the mean redshift of all galaxies ly- 
ing above 2 > 1.5 galaxies. This further demonstrates the 
low-2 peak is a low-significance contamination. A further 
examination of the photometric redshift estimation for such 
objects in the COSMOS field seems worthwhile in view of 
these results. 

5.3 Lensing strength dependence on magnitude 

As another check, we plot F, the py-amplitude ratio, vs. z'- 
band magnitude for A370 selected samples, to see if there is 
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Figure 5. Left: B] ~ Bq vs. Rq — z' CC diagram, showing the distribution of galaxies in COSMOS field. Applying the same CC-cuts as 
our selection for A370, we define galaxy samples: foreground (gray), orange (orange), green (green), red (red), blue (blue), and dropout 
(magenta) galaxies, in order to estimate the mean redshifts and mean depths of these samples using COSMOS photo-z's. Right: Average 
COSMOS redshift in CC bins. Overlaid are the boundaries of the foreground (thin dashed gray line), orange (thick dotted-dashcd orange 
line), green (thin solid green line), red (thick dashed red line), blue (thin dottod-dashed blue line) and dropout (thick solid magenta 
line) samples selected in the COSMOS field according to A370 CC-cuts. Evidently, different regions of CC space correspond to different 
redshift populations of galaxies. Most notably. The top left corner denotes dropout galaxies of 2 > 3.5, corresponding well to our selection 
of highly distorted dropout galaxies. 
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Figure 8. F, the qt amplitude ratio, vs. magnitude (z' band) for 
A370 orange, green, red, blue and dropout galaxy samples (top to 
bottom). Overlaid is the lensing depth, di^/da, vs. magnitude cal- 
culated from COSMOS photo-z's for each of the samples (dashed 
black line) with 1 ~ a confidence levels (shaded region). 



any trend of the lensing amplitude with magnitude. We also 
plot the mean di^/ds as a function of magnitude (all values 
normalized to the reference background values) calculated 
in independent magnitude bins for each sample using the 
COSMOS photo-z catalogue. This serves as a further con- 
sistency check. Interestingly, for the blue galaxies (and the 
dropout galaxies to some extent), the mean signal seems to 
drop slightly with fainter magnitudes. This trend is some- 
what counter-intuitive, since we expect that fainter galaxies 
will be at higher redshifts, and therefore have on average a 
higher signal. 

The diminished WL signal could possibly hint at a prob- 
lem with estimating the WL signal from faint blue galaxies, 
which are in general quite irregular in morphology, and em- 



pirical simulations with higher space based resolution can be 
made to examine this better. The declining trend of the pre- 
dicted rfis/rfs in the case of blue galaxies, based on the COS- 
MOS photo-z estimation, shown in Fig. 8 (dashed curves), 
could possibly point to a miss-classification of blue galaxies 
with photo-z methods, a well known problem for blue galax- 
ies, or even simply a limiting magnitude beyond which the 
photo-z method fails in this catalogue. We set our magnitude 
limits conservatively to minimize these possible problems, 
with 19 < z' < 23 for the orange sample, 20 < z' < 24.5 
for the green sample, 23 < z' < 25 for the blue sample, 
22 < z' < 26 for the red sample, and 24.5 < z' < 26.5 
for the dropout sample, in the case of A370. Similar exam- 
ination and limits were also applied for ZwC10024-|-17 and 
RXJ1347-11, with results discussed below. 



5.4 Lensing strength vs. redshift 

We may now finally combine the WL distortion information 
with the redshift information for all the samples of back- 
ground galaxies defined above. First, we plot the WL ampli- 
tude, a, as a function of redshift in Fig.9 for the three clus- 
ters - A370 (top), ZwC10024+17 (middle) and RXJ1347-11 
(bottom). A clear trend of increasing WL amplitude is seen 
with redshift, especially in the case of A370. A null result 
is easily excluded with low significance - /N — 91/8 for 
A370, x^/N = 34/4 for ZwC10024-|-17 and xV^ = 42/3 for 
RXJ1347-11, for a case of non-increasing horizontal line. In 
order to compare the resulting trend with the cosmological 
trend, we plot F, the jT-amplitude ratio, against median red- 
shift, (za) , in Fig. 10 for each of the clusters examined - A370 
(top), ZwC10024-|-17 (middle) and RXJ1347-11 (bottom). 
The horizontal bars show the width of the source redshift 
distribution, given as the 68% range of data about the me- 
dian redshift. Each point in a figure represents an indepen- 
dent sample, where the first point represents the foreground 



© 0000 RAS, MNRAS 000, 000-000 



10 Medezinski et al. 2010 



A3 70 



A370 





ZwC10024 



ZwC10024 




2 2.5 
<z> 



RXJ1347 




ACDM; X 'N=3.7/5 

EdS; x^'N=4.2/6 

n =0,ii =0; x^/N=3.2/5 



2 2.5 
<z> 

RXJ1347 



3.5 




ACDM; X /N=0.95/4 

EdS; x^/N=1/4 

n =0,£i =0; x^/N=0.8e/4 



1.5 



2.5 



Figure 9. gx amplitude, a, vs. rcdsiiift for A370 (top) brigiit 
(stars) and faint (circles) samples, ZwC10024+17 (middle) and 
RXJ1347-11 (bottom). A trend of higher WL amplitude with red- 
shift is seen. 



sample in front of each cluster, and the other points repre- 
sent background galaxy samples. For each sample, we also 
calculate the median lensing distance ratio, (dis/da), for each 
cosmological model - ACDM (empty circles), Einstein-de 
Sitter (crosses), and an empty universe (empty squares), us- 
ing the COSMOS photo-z measurements of galaxies within 
the same CC-magnitude boundaries defined for each back- 
ground sample and thus obtain the predicted depth. We 
interpolate between these discrete predicted values to pro- 



Figure 10. gT amplitude ratio, F, vs. redshift for A370 (top) 
bright (stars) and faint (circles) samples, ZwC10024-|-17 (mid- 
dle) and RXJ1347-11 (bottom). Also plotted is the lensing 
depth, dig/de, vs. redshift for different cosmologies - ACDM (cir- 
clesH-solid line), Einstein-de Sitter (crosses-|-dashed lino), and an 
empty Universe (squares-|-dashed-dotted line) estimated using the 
COSMOS photometric redshift catalogue. Horizontal bars repre- 
sent the width of the source redshift distribution, given as the 
68% range of data about the median redshift. 
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vide the theoretical relation, dis/ds{z) for each cosmological 
model. We can thus compare how the WL strength F agrees 
with predicted dis/ds. It is quite evident, especially in the 
case of A370, where we have the best dataset and therefore 
more data points to compare with, that the WL amplitude 
agrees well with the theoretical relations for dis/ds as a func- 
tion of redshift. For the other two clusters, ZwC10024-|-17 
and RXJ1347-11 the data is more shallow and hence noisier 
as can be seen in Fig. 10 (middle and bottom panels, re- 
spectively) but consistency is also found within the errors. 
Note that in the case of RXJ1347-11 the colour separation 
is not as good, since the bluest band available is V, rather 
than B, limiting the separation in depth of the background 
galaxy populations. 



6 CONSTRAINING COSMOLOGICAL 
PARAMETERS 

The clear detection here of the distance-redshift relation 
from our WL analysis of A370, prompts the question of 
how many such clusters would be required in order to pro- 
vide a useful cosmological constraint. This issue has been 
explored more generally in the context of planned field and 
cluster surveys by Jain & Taylor (2003); Taylor et al. (2007); 
Kitching et al. (2007). Taylor et al. (2007) present a detailed 
analysis of the sensitivity of cosmological cluster surveys to 
the ratio of shear values measured in independent redshift 
bins, finding that a large fraction of the potential integrated 
signal on the sky is contributed by abundant small, cluster 
mass range (M « lO^^M©), but that a large contribution 
also comes from the largest clusters, like those studied here. 

Here we use an order-of-magnitude calculation to esti- 
mate how well our newly approved MCT/CLASH* survey 
(P.L M. Postman) can do in this context, for which we aim 
to complete very high quality WL data for approximately 25 
massive clusters, similar in quality to the BjRcz' imaging 
of A370. 

As defined above in Eq. (12), F, is a shear ratio statistic 
between any two independent redshift bins summed over Nd 
clusters, given by: 



P . . ^ 7i ^ r[x{zj)]r[x{zi) - Xjzi)] 
~ 7i r[x{zi)]r[x{zj) - xizi)] 



(16) 



where r = r(x) is the comoving angular diameter distance, 
x(z) is the comoving distance and z; is the redshift of the 
lens, and F scales with the dark energy equation of state 



parameter, u;, as F 



(Taylor et al. 2007). The 



fractional error on w is given by (Taylor et al. 2007), 

a™ _ 2 /dlnF~ 
w 'yx \d\nw , 



(17) 



where 7t is the typical mean tangential shear of each cluster, 
and (Te = 0.3 is the measured intrinsic scatter in galaxy 
ellipticity per mode (KSB), and A*';, is the total number of 
galaxies summed up behind all the clusters co-added for this 
purpose. 

Assuming 7t ~ 0.05 and ~ 0.6 x 10®, summed over 
the available background for 25 clusters, (taking A370 as 



http : //www. stsci . edu/~postman/CLASH/ 



our guide to the number of background galaxies detected 
per cluster) we find from Eq. (17) the expected precision on 
w from our sample is Aw/w « 0.8. While this seems a rela- 
tively large uncertainty, the first application of the method 
by Kitching et al. (2007) suggests that the error distribution 
is non-Gaussian, with a rather sharp cut-off at high values 
of w placing a relatively tight upper limit. Other geometric 
probes currently do not succeed much better than this indi- 
vidually, e.g., from Baryon acoustic oscillations and SN-Ia, 
Aw « 0.3. Furthermore, the shear-ratio test has a different 
degeneracy with respect to the cosmological parameters to 
other probes, making even a crude measurement worthwhile. 

This constraint may be improved upon by conducting a 
more careful analysis incorporating a proper likelihood func- 
tion with reasonable priors for estimating the shear signal 
and redshift dependence. In addition, the photometry may 
be optimised for this purpose with emphasis on maximising 
dropout populations so that a wider redshift coverage may 
be achieved for better defining the distance-redshift relation. 



7 DISCUSSION & CONCLUSIONS 

Using deep observational data, the dependence of the ampli- 
tude of WL with source distance has been measured for in- 
dividual massive galaxy clusters using independent samples 
of foreground and background galaxies of differing depths, 
with a visible increasing trend. This is most clearly visible 
for A370 where we have Subaru imaging of relatively high 
quality m B, R and z' bands, allowing us to further subdi- 
vide these samples into independent bright and faint popula- 
tions. A general increasing trend with redshift is also seen for 
the other two clusters, albeit noisier than for A370. These 
clusters datasets are less deep, and with a reduced colour 
coverage (only V,R,z' in the case of RXJ1347-11), demon- 
strating the advantage of depth and the wide bandwidth 
coverage of A370. Small number statistics and possible di- 
lution, especially in the case of the dropout samples, may 
lead to large uncertainties and underestimated values. Our 
photometry comprises only three optical bands per cluster 
and so we do not rely on photo-z estimates, but instead we 
determine the depth of these background populations with 
reference to the very weU studied COSMOS and GOODS- 
MUSIC fields, by applying CC and magnitude cuts equal to 
that of our background populations. 

For A370, the trend of increasing WL-amplitude with 
redshift uncovered here follows the expected form of the 
lensing distance-redshift relation but with uncertainties 
presently too large to distinguish between cosmologies. The 
encouraging results from the clusters examined here, most 
notably for A370, merit further application of this approach 
to a larger sample of clusters. 

The recently approved MCT/CLASH program wiU ob- 
serve 25 clusters with FST ACS/WFC3, most of which have 
deep multi-colour Subaru imaging, so that we estimate a 
WL based precision on w of Aw « 0.6, but with a dif- 
ferent degeneracy relative to other probes, complementing 
existing methods. Combining this WL estimate with the 
distance-redshift relation from strong lensing will provide 
an enhanced geometric-based cosmological constraint. 

To extract the cosmological parameters from such accu- 
rate data will require further refinement of the method. We 
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must take account of the greater mean depth that lensing 
magnification generates whose effect on cluster lensing has 
been explored in some detail previously (Broadhurst et al. 
1995), and also a second order correction for the surface den- 
sity described in § 4.1. These effects, although small in terms 
of the lensing amplitude, are comparable with the relatively 
small differences of interest between competing cosmologies 
and thus must be explored in any serious study of cosmology 
with this method. 
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